extern doublereal pythag_(doublereal *a, doublereal *b);
